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^ ■ Inclusive production of the fo(980), f2(1270) and 0(1020) resonances has been studied in a 
P-i! sample of 4.3 million hadronic Z° decays from the OPAL experiment at LEP. A coupled channel 
I analysis has been used for the fo in simultaneous fits to the resonances in inclusive ir'^ir' and 
K"'"K~ mass spectra. Fragmentation functions are reported for the three states. Total inclusive 
rates are measured to be 0.141±0.007±0.011 fo, 0.155±0.011±0.018 fs and 0.091±0.002±0.003 
(j) mesons per hadronic decay. The production properties of the fo, including those in 
three-jet events, are compared with those of the f2 and 0, and with the Lund string model of 
hadron production. All measurements are consistent with the hypothesis that the fo(980) is a 
conventional qq scalar meson. 



Submitted to European Physical Journal C: Particles and Fields 



The OPAL Collaboration 



K. AckerstafF, G. Alexander J.AUison^^, N. Altekamp^, K.J. Anderson^, S. Anderson^^, 
S.Arcelli^, S.Asai^^, S.F.Ashby^, D.Axen^^, G. Azuelos^^-", A.H.Balli^ E. Barberio^ 
R.J.Barlow^^ R. Bartoldus^ J.R. Batley^ S.Baumann^, J. Bechtluft^^ T. Bchnke^ 
K.W.BelP, G.Bella^^^ S. Bentvelsen^ S.Bethkel^ S.Bettsl^ O. Biebel^^ A. Biguzzi^ 
S.D.Bird^^ V.BlobeF, I.J. Bloodworthi, M. Bobinski^o, R Bock^^ D. Bonacorsi^, 
M. Boutemeur^^, S.Braibant^, L. Brigliadori^, R.M. Brown^", H.J. Burckhart^, G. Burgard^, 
R. Biirgin^", P.Gapiluppi^, R.K. Carnegie^, A.A. Garter^^, J. R. Garter^, G.Y. Ghang^^, 
D.G.Gharltonl'^ D.Ghrisman^ RE.L. Glarkc^^, I. Gohcn^^, J.E. Gonboy^^, O.G.Gooke^ 
G. Gouyoumtzclisi'\ R.L. Goxc^ M.Guffiani^ S.Dado22, G. Dallapiccolai^, G.M. Dallavalle^, 
R.Davis^o, S. De Jong^^ L.A.del PozoS A. de Roeck^ K. Desch^ B. Dicncs'^'^'^ M.S. Dixit^ 
M.Doucet^^ E.DuchovnP, G.Duckeck^^, I.R Duerdoth^^ D.Eatoughl^ P.G. Estabrooks^ 
E. Etzion23, H.G. Evans^ M. Evans^^ F. Fabbri^, A. Fanfani^, M. Fanti^ A.A. Faust^o, L. Feld^, 

F. Fiedler^^, M. Fierro^, H.M. Fischer^ I. Flrf, R. Folman^^, D.G.Fong^^ M. Foucher^^ 

A. Furtjes^ D.I. Futyan^^, P. Gagnon^ J.W. Gary^ J. Gasconl^ S.M. Gascon-Shotkin^^ 
N.I.Geddes^o, G. Geich-GimbeI^ T. Geralis^o, G. Giacomelli^ P. GiacomeIIi^ R. Giacomelli^, 
V.Gibson^ W.R. Gibson^^^ D.M. Gingrich^O'", D. Glenzinski^, J. Goldberg^^, M. J. Goodrf, 

W.Gorn^ G. Grandi^, E.Gross^^ J. Grunhaus^^, M. Gruwe^^, G.Hajdu^^, G.G. Hanson^^ 
M. Hansroul^, M. Hapke^^, G.K. Hargrove'^, P.A.Hart^, G. Hartmami'^ M. Hauschild^ , 

G. M. Hawkes^, R. Hawkings^^, R.J. Hemingway^, M. Herndon^'', G.Hcrtcn^°, R.D.Heuer®, 

M.D. Hildreth^ J.G. HiII^ S.J. Hillieri, P.R. Hobson^^, A. Hocker^ R.J.Homer^ 
A.K. Honma^^'", D. Horvath^^-^, K.R. Hossain^", R. Howard^^ P. Hiintemeyer^^, 
D.E.Hutchcroft^ P. Igo-Kemenes^\ D.G.Imrie^^ K.Ishii^^, A. Jawahery^"^, P.W. Jeffreys^°, 
H. Jcrcmie^^, M. Jimack^, A. Joly^^, G.R. Joncs^, M. Jones^, U. Jost^^, P. Jovanovic^, 
T.R. Junk^, J.Kanzaki^^, D.Karlen^, V. Kartvelishvili^^, K. Kawagoe^"^, T. Kawamoto^^, 
P.I.KayaF, R.K. Keeler^^, R.G. Kellogg^^ B.W. Kennedy^o, J.Kirk^^, A.KIier^e, S.KIuth^, 
T.Kobayashi^^, M.KobeP°, D.S.Koetke^ T.P.Kokott^ M.Kolrep^^, S. Komamiya^^, 
R.V.KowalewskP, T.Kress^^ P.Krieger^ J. von Krogh^\ P.Kyberd^^, G.D. Lafferty^^, 
R. Lahmann^^, W.P. Lai^^, D.Lanskc^^, J. Lauber^^, S.R. Lautenschlager^^, I.Lawson^^, 
J.G.Laytcr^ D.Lazic^^, A.M.Lec^i, E. Lefebvrel^ D.Lellouch^^, J. Lettsl^ L. Levinson^^, 

B. List^ S.L.LIoydi^^ F.K. Loebingeri^, G.D.Long^^, M.J.Losty^ J. Ludwig^o, D.Lui^^ 
A.Macchiolo^, A. Macpherson^°, M. Mannelli^, S. Marcellini^, G. Markopoulos^^, G.Markus^, 

A.J. Martin^^, J.P. Martin^^, G. Martinez^'^, T. Mashimo^^, P. Mattig^^, W.J. McDonald^°, 
J.McKenna^s, E.A. Mckigney^^ T. J. McMahon\ R.A. McPherson^^ F.Meijers^ S.Menke^, 
F.S.Merritt^ H. Mes^ J. Meyer^^, A. Michelini^ S. Mihara^^, G. Mikenberg^^, D.J. Miller 
A. Mincer^^'^, R. Mir^^, W. Mohr^°, A. Montanari^, T.MorP, S. Miliar a^^ K. Nagai^^, 
LNakamura^^, H.A.NeaF^ B.Nellen^ R.Nisius^ S.W. 0'Neale^ F.G. Oakham^ F.Odorici^, 
H.O.Ogren^^^ A.Oh^^, N. J. 01dershaw^^ M.J. Orcglia^ S.Orito^^, J. Palinkas'^'^ '^, 
G.Pasztor32, J.R. Pater^^ G.N. Patrick^o, J.Patt^o, R. Perez-Oclloa^ S. Petzold^^ 
P. Pfeifenscllneiderl^ J.E. Pilclier^ J. Pinfold-^o, D.E.Plane^ P. Poffenberger^^ , B.Poli^, 
A. Posthaus^, G.Rembser^, S. Robertson^^ S.A. Robins^^ N.Rodning^o, J.M.Roney^^, 
A.Rooke^^, A.M. Rossi^, P. Routenburg^^, Y. Rozen^^, K.Runge^°, O. Runolfsson^, 
U.RuppeP^ D.R.Rust^^^ K.Sachs^", T. Saeki^^ O.Sallr=^^ W.M.Sang^^, E.K.G. Sarkisyan^^, 
G.Sbarra^s, A.D. Schaile^^ 0. Scliaile^^, F. Scharf^, P. Scharff-Hansen^, J. Schieck^^ 
P. Schleper^^ , B. Schmitt^, S.Schmitt^^, A. Schoning^, M. Schroder^, M. Schumacher^, 
G.Schwick^ W.G.ScotP, T.G.Shears^ B.G.Shen^ G.H. Shepherd-Themistocleous^ 



1 



P.Sherwood^^ G.P.Siroli^, A. Sittler^^, A. Skillmanis, A. Skuja^^ A.M.Smith^ G.A. Snow^^ 
R. Sobie^^, S. Soldner-Rembold^'', R.W. Springer^'^, M. Sproston^°, K. Stephens^^, J. Steuerer^^, 
B. Stockhauscn'\ K. Stoll^'', D. Strom^^, R. Strohmcr'^^, P. Szymanski^'^ , R. Tafirout^^, 
S.D. Talbot^, P. Taras^^, S. Tarem^^, R. Teuscher^, M. Thiergen^°, M. A. Thomson^, E. von 
Torne^, E.Torrence^, S. Towers^, I.Trigger^^, Z. Trocsanyi^^, E. Tsur^^, A.S. Turcot^, 
M.F.Turner- Watson^, I.Ueda^^, P.Utzat", R.Van Kooten^^ P. Vannerem^°, M. Verzocchi^°, 
P.Vikasl^ E.H. Vokurkal^ H.Voss^, F. Wackerle^o, A.Wagner^^, C.P.Ward^ D.R.Ward^ 
P.M.Watkins\ A.T.Watson\ N.K.WatsonS P.S.Wells^ N.Wermes^, J.S. Whitens, 
G.W. Wilson^^, J.A.Wilson^, T.R.Wyattl^ S. Yamashita2^ G.YekutielP, V. Zacek^^^ 

D. Zer-Zion^ 



^School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
^Dipartimento di Fisica dell' Univcrsita di Bologna and INFN. 1-40126 Bologna, Italy 
^Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
^Department of Physics, University of California, Riverside CA 92521, USA 
^Cavendish Laboratory, Cambridge CBS OHE, UK 

^Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 
Ontario KIS 5B6, Canada 

''Centre for Research in Particle Physics, Carleton University, Ottawa, Ontario KIS 5B6, 

Canada 

^CERN, European Organisation for Particle Physics, CH-1211 Geneva 23, Switzerland 
^Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 
USA 

^^Fakultat fiir Physik, Albert Ludwigs Universitat, D- 79104 Freiburg, Germany 

^^Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 

^^Indiana University, Department of Physics, Swain Hall West 117, Bloomington IN 47405, 

USA 

^^Queen Mary and Westfield College, University of London, London El 4NS, UK 
^^Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D-52056 

Aachen, Germany 

^^University College London, London WCIE 6BT, UK 

^^Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 
^^Department of Physics, University of Maryland, College Park, MD 20742, USA 
^^Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 

Canada 

^^University of Oregon, Department of Physics, Eugene OR 97403, USA 
^"Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
^^Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 
^^Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
^^International Centre for Elementary Particle Physics and Department of Physics, University 
of Tokyo, Tokyo 113, and Kobe University, Kobe 657, Japan 

^^Institute of Physical and Environmental Sciences, Brunei University, Uxbridge, Middlesex 
UB8 3PH, UK 

^^Particle Physics Department, Wcizmann Institute of Science, Rehovot 76100, Israel 
^^Universitat Hamburg/DESY, II Institut fiir Experimental Physik, Notkestrasse 85, D-22607 
Hamburg, Germany 



2 



^^University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, Canada 
^^University of British Columbia, Department of Physics, Vancouver BC V6T IZl, Canada 
^"University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
3iDuke University Dept of Physics, Durham, NC 27708-0305, USA 

^^Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hungary 
^^Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 
^^Ludwigs-Maximilians-Universitat Miinchen, Sektion Physik, Am Coulombwall 1, D-85748 
G arching, Germany 



" and at TRIUMF, Vancouver, Canada V6T 2A3 
* and Royal Society University Research Fellow 
and Institute of Nuclear Research, Debrecen, Hungary 

and Department of Experimental Physics, Lajos Kossuth University, Debrecen, Hungary 
^ and Department of Physics, New York University, NY 1003, USA 



3 



1 Introduction 



Inclusive production of mesons and baryons in TP decay has been studied extensively and 
the results have provided valuable input to the theory and phenomenology of parton hadroniza- 
tion. Production of the fo(980) is of particular interest because although it is well established 
experimentally as a scalar (J^*-^ = O"*"^) state, its precise nature has long been uncertain. In- 
deed, there is still vigorous debate about the identity of the mesons which comprise the 
lowest lying scalar nonet of flavour SU(3). Compared to the expectations for a conventional 
meson at a mass of around 1 GeV, the fo has a markedly small total width, a relatively large 
coupling to KK, and a partial width to 77 which is an order of magnitude below theoretical 
expectation. 

A number of suggestions have been made as to the nature of the fo. Jaffe and Johnson ||^ 
performed a bag model calculation to suggest that it could be a "cryptoexotic" qqqq state. 
Weinstein and Isgur 0, using a potential model of qqqq states, showed that the fo could be 
explained as a loosely bound KK system, a so-called KK molecule. Gribov [^,|^ has proposed a 
theory of confinement in QCD, in which the fo plays the role of a novel "vacuum scalar" state, 
a bound state of a quark and antiquark with negative kinetic energy, interacting repulsively 
to give a state of positive total energy. Ishida et al. have proposed an interpretation as 
a hybrid meson with a massive constituent gluon, while a scalar glueball has been suggested 
by Robson 0. An analysis by Close and Amsler of the scalar states suggests that the 



most likely candidate for the lowest-lying scalar glueball is the fo(1500); in this model, the 
fo(980) is a left-over state which cannot be accommodated in the scalar nonet. Lattice QCD 



calculations |Tl[] also find the lightest scalar glueball mass to be around 1.5 GeV so that the 
glueball interpretation of the fo is now out of favour. Tornqvist [|T^], using a unitarized quark 
model to analyse all of the identified states in the O"*""*" sector, has concluded however that the 
fo can fit in as a conventional qq meson in the lowest 0"'""'" multiplet. Similarly, an analysis by 
Zou and Bugg [|13| of all available high-statistics vr+vr^ and KK scattering data concluded that 



the fo could be interpreted as a meson, and recent work by Anisovich and Sarantsev |TJ] has 
come to the same conclusion. 

One aim of the present study is to look for features of fo production in TP decay which 
may help to elucidate its nature. The vacuum scalar states of Gribov's theory are expected 
to be compact objects with distinctive production properties, and suggestions have been made 
by Close et al. for a number of experimental tests. In Z° decay, the signature would be a 
relatively larger yield in low multiplicity events and in events where the fo is isolated in rapidity. 
If, on the other hand, the fo has a significant gluonic content, its production could be enhanced 
in gluon jets. In contrast, if it is principally a conventional meson, its production properties 
may be unremarkable when compared to those of other similar states. There are presently no 
predictions for production rates of qqqq states or KK molecules in Z° decay. 

Two approaches are taken in the present analysis: a comparison of the features of fo produc- 
tion with those of two neutral isoscalar mesons, and a comparison with the JETSET 7.4 



implementation of the Lund string model of hadronization [|T6|, within which the fo is treated 
as a conventional meson. The features of fo production are measured along with those of the 
f2(1270) and the 0(1020) mesons. The latter is close in mass to the fo, while the former is an 
established meson with the qq in a state. In the conventional meson interpretation, the fo 
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would be a ^Pq state. Since the mass difference is relatively small, the fg and £2 may then be 
expected to have similar production properties. The Lund string model, as implemented in the 
JETSET 7.4 Monte Carlo program, is highly successful in describing many features of hadronic 
Z° decay events. These features include global event properties related to the perturbative 
phase of the initial parton shower, as well as details of the production of many different species 
of hadrons during the nonperturbative hadronization phase. In short, the model provides a 
well-tested and reliable picture of particle production in hadronic T? decays |]T^. Within the 



JETSET model, the fo is treated as a scalar meson composed of uu and dd pairs, and the 
relative production rates of ^Pq and ^P2 mesons are determined by a variable parameter. A 
comparison of the data and the results of JETSET could therefore provide information on the 
nature of the fo. 

The only reported measurements of the fo and f2 at LEP were made by DELPHI as part 
of a general study of light resonances [l^ using relatively low statistics. The fragmentation 



functions and total rates were reported over restricted ranges of momentum. In a similar study 
of resonance production by ALEPH [l^, the fo was included in fits to tt+tt" mass spectra, 
but no rates were reported. Because the fo has a significant coupling to KK and peaks below 
threshold for this channel, its 7r"^7r~ mass spectrum will not in general exhibit the simple Breit- 
Wigner shape assumed in these two previous analyses. In the present study, a coupled-channel 
analysis is done, including tt+tt^ and K"'"K~ data, in which proper account is taken of the 
opening of the KK channel in fo decay. 



2 The OPAL detector and data samples 



The OPAL detector is described in |]20|. For the present analysis, the most important com- 



ponents were the central tracking chambers which consist of two layers of silicon microvertex 



detectors ||21|, a high-precision vertex drift chamber, a large- volume jet chamber, and a set of 
drift chambers (the 2;-chambers) which measure the coordinates of tracks along the direction 
of the beam. The OPAL coordinate system is defined with the ^-axis following the electron 
beam direction; the polar angle Q is defined relative to this axis, and r and are the usual 
cylindrical polar coordinates. The central chambers lie within a homogeneous axial magnetic 
field of 0.435 T. Charged particle tracking is possible over the range | cos6'| < 0.98 for the full 
range of azimuthal angles. The OPAL jet chamber is capable of measuring specific energy loss, 
d£^/da;, with a resolution, cr(d£'/dx)/(di?/da;), of 3.5% for well-reconstructed, high-momentum 
tracks in hadronic events [p2| . 

The present analysis used the full OPAL sample of 4.3 million hadronic T? decays recorded 
at LEP 1 between 1990 and 1995. To correct for losses due to the acceptance and efficiency of 
the experiment and the selection procedures, and also to provide signal and background shapes 
for fits to the data mass spectra, 6 million Monte Carlo events were used, which had been 
generated using JETSET 7.4 and processed through a full simulation of the experiment 



and the data reconstruction and analysis. This 'detector-level' Monte Carlo sample was also 
used for comparison with various features of the experimental data. The JETSET version was 
tuned |4| using OPAL data on event shape distributions, fragmentation functions of vr^, K^, 



p/p and A, and LEP data on total inclusive multiplicities for 26 identified hadron species. 
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A detailed description of the selection of hadronic Z° decay events in OPAL is given in |25 . 
For the present analysis, tracks in the selected events were required to have: a minimum 
momentum transverse to the beam direction of 150 MeV/c; a maximum momentum of 1.07 x 
-E'beam, based on the momentum resolution of the detector; a distance of closest approach to 
the interaction point less than 5 cm in the plane orthogonal to the beam direction, and the 
corresponding distance along the beam direction less than 40 cm; a first measured point within 
a radius of 75 cm from the vertex; and at least 20 hits available for measurement of specific 
energy loss, dE/dx. 

Kaons and pions were identified using the dE/dx measurements. For each track, a 
probability (weight) was formed for each of the stable particle hypotheses (e, fi, vr, K and p). A 
track was identified as a pion or a kaon if the appropriate weight was above 5% and was larger 
than the weight for each of the other stable hadron hypotheses. Between momenta of 0.8 and 
2.0 GeV/c, the vr, K and p bands overlap in dE/dx, leading to considerable ambiguity among 
hypotheses. Therefore, no tracks were identified as kaons in this momentum range although, 
since most particles are pions, pion identification was still allowed. 

With the event and track selection cuts previously described, fo and f2 decaying via tt+tt" 
were identified with an efficiency of around 40% over the whole momentum range. The — *■ 
K"'"K~ efficiency typically varied between 15 and 20%, although in the (j) momentum range from 
1.4 to 3.6 GeV/c it fell below 5% because of the dE/dx cross-over region. The mass resolution 
at 1 GeV varied with momentum, from 15 to 20 MeV for tt+tt" and from 2.5 to 4.5 MeV for 
K+K". 



3 Data analysis 

3.1 Inclusive two-particle mass spectra 

For simultaneous fits to tt+tt^ and K'^'K" mass spectra in both the real and simulated data 
samples, spectra were formed for inclusive n'^n^ , tt^tt^ and K+K" systems in bins of each of 
the following variables: 

• scaled momentum, Xp{= pZ-Ebcam), of the two-particle systems (nine bins) 

• rapidity, y, of the two-particle systems with respect to the event thrust axis (six bins) 

• rapidity gap. Ay, between the two-particle systems and the closest single charged particle 
(six bins) 

• multiplicity, rich, of charged tracks in the event (six bins) 

The last three of these variables were chosen specifically for tests of the nature of the fo, as 
discussed in section 1. To account for a large part of the combinatorial backgrounds in tt+tt", 
mass spectra were formed by subtracting the tt+tt^ spectra from those for 7r~^n~ . Because the 
combinatorial backgrounds in K+K" were relatively small, and Bose-Einstein correlations could 
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affect the K^K^ spectra near threshold, no subtraction was done for the K+K~ mass spectra. 
For the Monte Carlo sample, separate spectra were also made for the most important states 
contributing to the mass spectra, using information on the origin of each track at the generator 
level. 



3.2 Selection of three- jet events 

To investigate possible differences between production in quark and gluon jets, the Durham jet 



finder |23] was used to identify a sample of three-jet events. The cut-off value ?/cut was set to 
0.005, and the jet-finding was done using charged tracks. For each candidate three-jet event, 
the angle between the two lowest-energy jets was required to be greater than 30°, and in order 
to ensure well-reconstructed, planar events, the sum of the interjet angles was required to be 
larger than 358°. The jet energies were then reconstructed using the interjet angles, assuming 
massless kinematics. Each jet was required to contain at least two charged particles and more 
than 5 GeV of energy. With these cuts, 24% of all events were selected as three-jet events. 
Monte Carlo studies have shown that in such events, 95% of the highest-energy jets are due to 
quarks, while the lowest-energy jet is a gluon jet with an approximately 80% probability. Mass 
spectra were formed, using only tracks assigned to the same jet, in three bins of E/Ejet where 
E is the energy of the two-particle system and E'jet is the energy of the jet. 



3.3 Fit procedures 

For each bin of the above kinematic variables, the 7r"'"7r~ and K+K~ mass spectra were fitted 
simultaneously, using a minimum fit, to a sum of contributions given by: 

fifn^+n") = aj^"" |Afo(m^^)p + af^BWf^ + al'^J B^+^- (1) 

/(mx+K-) = ^"|A,(^kk)|' + a^BW,^ + ag^/'^K+K- (2) 

In equations (|l]) and (|), the a terms represent the intensities to be fitted, A^^ is the amplitude 
for fo(980), BW are Breit-Wigner functions, and the B represent background functions. The 
are related by 

af = 0.75^^^aF;^" 

where eK+K-/e7r+7r- is the ratio of the efficiency to reconstruct fo in K+K" relative to tt+tt" 
(which varies with the bin of the kinematic variable), and the factor 0.75 comes from Clebsch- 
Gordan coefficients. Following Flatte |2^, the coupled-channel amplitudes for fo(980) decay 
via TTTT and KK were taken to be: 



^fol^l^TrJ- 2_^2 _irn (V ^V...A ^' 
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KiniKK) = 2 ^2 (r ) 

^0 ~ ^KK ~ U TTTT + J- KKJ 

Here rriQ is the resonance mass, and the partial widths F are related to the coupling constants 
g via 

T^^ = g^\i—^ -ml and Tkk = fi-Ky - "^k 

The fact that Fkk is imaginary below KK threshold leads to distortion of the vrvr mass spectrum 
from a simple Breit-Wigner shape. Interference between fg and 7r"'"7r~ backgrounds in Z° decay 
should be negligible and was not included in the fits. In determining the meson rates, the 
possible presence of other f-meson resonances coupling to tt+tt^ has been neglected. 

The experimental mass resolution was treated differently for the three states. The fo inten- 
sities derived from equation (3) were folded with the detector resolution separately for each bin 
of Xp, Ay, ricYi and E/Ejf,^. For the broad f2, a relativistic D-wave Breit-Wigner was used, 
with mass and width fixed to the Particle Data Group (PDG) values 0. For the (p, where 
experimental mass resolution is particularly important, the shapes for the fits were taken from 
the detector-level output of the simulation, accounting automatically for variations of the mass 
resolution. The background shapes for the fits were also taken from the detector-level Monte 
Carlo: for tt+tt", all pairs identified as tt+tt^ were taken, except those from the fo and f2; and 
for the K+K~, all identified combinations were used except those due to fo and 0. In the fits, 
the background shapes were given additional freedom, with the tt~^tt~ contribution multiplied 
by a term (m"+^_ + a2'^"+^-) and the K+K" shape by (1 -|- a^rriK+K-), with the a being 
variable parameters. The p° resonance in the JETSET simulation is truncated at ±2.5F; to 
allow for this, the p*^ intensity above 1.15 GeV was taken as a straight-line extrapolation from 
the contribution at 1.15 GeV falling to zero at 1.5 GeV. 

The mass ranges of the fits were 0.82 to 1.5 GeV for tt+tt^, and from threshold to 1.18 GeV 
for K+K", these limits being chosen mainly to avoid the peak region in tt+tt" and the region 
of its reflection in K+K" (caused by particle misidentification). Figure shows as an example 
the sum of the fitted mass spectra over the range Xp > 0.14, indicating clearly the signals due 
to the three resonances. 

In fits to the total mass spectra for Xp > 0.14 with the fo parameters allowed to vary, the mass 
was well-constrained while the couplings g were poorly determined: ttiq = 0.957 ± 0.006 GeV, 
S^TT = 0.09 ± 0.40 and gK = 0.97 ± 0.82. This is essentially due to the limited mass resolution 
and the large backgrounds in the tt+tt^ channel; the K+K" data are relatively insensitive to 
the fo parameters. The fitted mass value is in good agreement with that obtained in the high- 



statistics analysis made by Zou and Bugg [|l3l of tttt — KK elastic scattering phase shifts and 



TTTT and KK mass spectra from central production in pp collisions. The mo and g values were 
therefore fixed to the Zou and Bugg values, mo = 0.9535 GeV, (?7r=0.111 and gK = 0.423. These 
couplings correspond to a branching ratio, BR(fo vrvr) = 0.80, in agreement with the PDG 
value. Nevertheless, there is much experimental and theoretical uncertainty in the resonance 
parameters appropriate for the fo. Indeed the PDG quotes the total width as 40 to 100 MeV. 
The results should therefore be taken as model-dependent measurements, which assume the 
Flatte parametrization with the parameter values of Zou and Bugg, and no interference with 
background. Fits to the 7r"'"7r~ and K"''K^ spectra are shown in figures § and ^ for six bins of 
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Xp] it can be seen from the figures that the parametrization works well in the decay data 
over the entire kinematic range. 



3.4 Fragmentation functions and total rates 

The rates obtained from the fits in bins of Xp are given in table |I|. Figure ^ shows these 
measured fragmentation functions, (l/(7h)do"/dXp (where cxh is the total hadronic cross section), 
along with curves from the JETSET 7.4 Monte Carlo generator, normalized for each particle 
to the measured total rate seen in the data. Only statistical errors are shown for the fo and 
f2; the dominant systematic errors are correlated over the Xp bins and are discussed below in 
section |3.5| . There is clearly good agreement for all three particles between the shapes of the 
momentum distributions in the Monte Carlo model and those in the data. 

The total inclusive rates have been measured, by integrating the fragmentation functions, 
to be: 

0.141 ±0.007 ±0.011 fo(980) 
0.155 ±0.011 ±0.018 f2(1270) 
0.091 ±0.002 ±0.003 0(1020) 

per hadronic Z° decay. The first errors are statistical and the second systematic (discussed 
below). These rates compare with values of 0.154 ± 0.025 fo and 0.240 ± 0.061 f2 mesons. 



Xp. The measurement for the 
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28 



results, but is three standard deviations below an ALEPH |2^ measurement. The present 



measurement has the smallest statistical and systematic errors. 



3.5 Systematic errors 

Table |^ gives a summary of systematic errors on the meson rates. To estimate systematic 
errors for the fo, the uncertainties in the resonance mass and the couplings were considered. 
The value of mo was varied from 0.951 to 0.963 GeV (the range obtained when the mass was 
allowed to be free in the fits) while keeping the coupling constants fixed. The maximum total 
variation, 2.5%, was taken as the systematic error. Next uiq and gK were fixed and the coupling 
(^TT was allowed to vary freely in the fit. This gave a 7% change in the total rate, which was 
conservatively taken as the systematic error from this source. (If, instead, the errors in the tttt 
branching ratio from the PDG are used, the systematic error would be 3%.) Finally, mo and 
Qt^ were fixed and was allowed to vary. The fo rate changed by much less than the statistical 
error, while the rate increased by 1.5%; this was assigned as a systematic error on the 
measurement. Systematic errors from uncertainties in Monte Carlo modelling of the track cuts 
contributes 1.4% ||2^ for resonances decaying to two charged particles. The overall systematic 
error on the fo measurements was therefore 7.6%. 

For the f2 measurements, the main sources of systematic uncertainty also come from the 
resonance parametrization, although in this case it is safe to assume a normal Breit-Wigner 
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resonance. To take account of possible long tails in the line shape beyond the upper limit of the 
fits to the mass spectra, the relativistic D-wave Breit-Wigner used for the fits was integrated out 
to 2.2 GeV, five full widths above the nominal peak position. This resulted in an increase of 30% 
in the total intensity over the fitted value. Since the shape of such a resonance is uncertain so 
far from the pole position, one half of this extra contribution was added to each measurement of 
the differential cross section (the results in table |l| include these corrections), and a systematic 
error of 30%/v^T2 was assigned. The fits were repeated, varying the f2 mass and width by 
one standard deviation (using the PDG values) above and below their nominal values. The 
maximum change in the measured rate, 6%, was assigned as a systematic error. To account 
for the high mass tail of the p° above 1.15 GeV, the fits were repeated assuming two extreme 
possibilities: the default JETSET simulation, with the shape truncated at 1.15 GeV; and a 
constant contribution set to the level in the simulation just below 1.15 GeV. The maximum 
change, 3.2%, was taken as the systematic error. With an error of 3.1% on the branching ratio 
to n~^7r~ and 1.4% from the modelling of the track cuts, the overall systematic error on the f2 
measurements was therefore 11.6%. 

Systematic errors due to the simulation of the energy loss were measured by varying the 
assumed mean values of the theoretical dE/dx distributions for a given particle hypothesis, 
and the assumed resolution on the energy loss measurements. Studies of well-identified pions 
from Kg decays, protons from A decays and kaons from D° decays were used to place limits 
on the maximum possible deviations of these quantities, and the analysis was repeated, with 
the dE/dx weights of the tracks being recalculated each time. Systematic errors were assigned 
as the maximum measured deviations from the standard fit values. The resulting errors were 
negligible for the fo and f2, and contributed a 1.3% error on the total rate. Since these errors 
varied with momentum, they are also given in table |^ as uncorrelated systematic errors. Because 
the (f) is close to threshold, there is a small background due to conversion e"'"e~ pairs, particularly 
at low Xp (as seen in figure 3). In the standard fit, this was included in the single background 
term, i?K+K- of equation (§). When this component was allowed to vary independently, the 
total (p rate changed by 2%, which was taken to be the systematic error from this effect. As 
in ||2^, an error of 1% comes from uncertainty in the mass resolution and 1.4% from modelling 
of the track cuts. Finally an error of 1.2% was included for the (j) —>■ K+K^ branching ratio, 
bringing the overall systematic error to 3.5% (including the contribution from variation of Qk 
discussed above). 

In all of the measurements, the possible presence of other f-meson resonances (for example 
the fo(1360)) in the tt+tt^ mass spectrum has been neglected. While there is no evidence in the 
data for such states, and therefore no need to include them in the fits, their presence could in 
principle give rise to additional systematic errors, particularly for the f2. 



3.6 Production as functions of event multiplicity and rapidity gap 

As discussed in section |l|, the Gribov vacuum scalar states are expected to be produced prefer- 
entially in low-multiplicity events, and when isolated in rapidity relative to the other particles. 
Figure ^ shows the fits to the 7r"'"7r~ mass spectra in the fg mass region for the six bins of 
event charged-track multiplicity. Figure ^ gives, for each of the three resonances, the ratios 
of production in data relative to the detector-level JETSET sample (in which the resonances 
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are treated as conventional mesons) as a function of the multiplicity. Figure |^ shows the same 
ratios as a function of the magnitude of the rapidity difference between the resonances and the 
nearest charged particle. For each resonance the ratios have been normalized such that the 
weighted averages are unity. No evidence is seen for anomalously large production of fo(980) 
either at low multiplicity or at large rapidity gap, with the Monte Carlo model giving a good 
description. The data therefore do not provide any evidence to support the hypothesis that the 
fo is a manifestation of a vacuum scalar state. 

3.7 Production in quark and gluon jets 

The total rates in each of the jets in the selected three-jet events were measured by fitting to the 
mass spectra in the three bins of scaled energy, and summing the contributions. The jet-finding 
procedures described in section pl2| were also applied to the detector-level Monte Carlo sample, 
and the total number of each of fo, f2 and in each jet was counted. Figure ^ shows the ratios 
of total rates measured in the data relative to the detector-level Monte Carlo model for each 
of the three jet classes. The figure shows no significant evidence for any differences between 
production in the quark-enriched (high-energy) and the gluon-enriched (low-energy) jets. 



4 Summary and conclusions 

Fragmentation functions and total inclusive rates in Z° decay have been measured for three 
resonances, the fo(980), f2(1270)and (/>(1020), using the full LEP 1 statistics of OPAL. The 
fo and f2 measurements will provide input to understanding the physics of inclusive particle 
production in the P-wave meson sector, particularly in the context of the Monte Carlo models. 

The production characteristics of the fo show no significant differences from those of the f2 
and (j) mesons. In particular, the shapes of the fragmentation functions are similar for all three, 
and are well reproduced by the JETSET 7.4 Monte Carlo model, within which the fo is treated 
as a conventional scalar meson. The total inclusive rate for the fo is 0.141 ± 0.013 per hadronic 
Z° decay, similar to the value, 0.155 ± 0.021, measured for the f2. The fo rate is significantly 
larger than for the 0, 0.091 ± 0.004 per hadronic Z° decay. These features are consistent with 
the fo being, like the f2, a qq meson in the '^P state composed mainly of uu and dd. 

For all three states, the production rates relative to the Monte Carlo model have been 
measured as functions of the charged particle multiplicity of the event and the gap in rapidity 
to the nearest charged particle. The distributions are found to be flat in all cases. In particular, 
no evidence has been found for enhanced fo production at low multiplicities or at large rapidity 
gap. There is thus no evidence to identify the fo with the vacuum scalar state proposed by 
Gribov. 

Production in energy-ordered jets in three-jet events has been measured with a view to 
seeking differences between quark and gluon induced jets. No significant differences are seen 
between the data and the Monte Carlo model in the relative production rates. There is therefore 
no evidence for any enhanced gluon content in the fo. 
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In summary, all measured characteristics of fo(980) production in the Z° decay data of OPAL 

are consistent with its interpretation as a conventional scalar meson. Quantitative theoretical or 
phenomenological predictions for production of the types of states discussed in the introduction 
could enable more definite conclusions to be drawn from the Z° data. 
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Xp range 


fo(980) 


f2(1270) 


0(1020) 


0.00-0.06 


1.04±0.09 


1.00±0.14 


0.464±0.011±0.005 


0.06-0.12 


0.57±0.05 


0.69±0.08 


0.316±0.021±0.007 


0.12-0.14 


0.30±0.06 


0.41±0.09 


0.285±0.020±0.009 


0.14-0.16 


0.20±0.05 


0.25±0.08 


0.197±0.019±0.006 


0.16-0.20 


0.21±0.03 


0.27±0.04 


0.167±0.017±0.002 


0.20-0.25 


0.13±0.02 


0.22±0.03 


0.133±0.007±0.002 


0.25-0.35 


0.085±0.011 


0.091±0.016 


0.096±0.004±0.001 


0.35-0.50 


0.046±0.005 


0.035±0.008 


0.045±0.002±0.001 


0.50-1.00 


0.0079±0.0009 


0.008±0.001 


0.010±0.001±0.000 


0.00-1.00 


0.141±0.007 


0.155±0.011 


0.091±0.002±0.001 



Table 1: Measured differential cross sections, {1/ ah)da /dxp, and total inclusive rates for fo, 
f2 and (f) production. There are additional overall systematic errors of 7.6% for the fo, 11.6% 

for the f2 and 3.5% for the (f). The uncorrclated systematic errors for the are given in the 
table; the error for the Xp range 0.50 — 1.00 is smaller than 0.0005. The uncorrelated systematic 
errors for the fo and f2 are negligible in comparison with the statistical errors and the correlated 
systematic errors. 





fo(980) 


f2(1270) 


0(1020) 


Variation of mo 


2.5% 






Variation of 51-^ 


7.0% 






Variation of (^k 






1.5% 


f2 resonance line shape 




8.7% 




f2 mass and width 




6.0% 




BR(f2 ^ TT+TT-) 




3.1% 




p° resonance line shape 




3.2% 




<\E/<\x parametrization 






1.3% 


BR(0 ^ K+K-) 






1.2% 


K+K~ mass resolution 






1.0% 


Photon conversions 






2.0% 


Track cuts 


1.4% 


1.4% 


1.4% 


Total 


7.6% 


11.6% 


3.5% 



Table 2: Sources of systematic errors on the total inclusive meson rates. Errors on the rates 
from uncertainties in the parametrization of dE/dx are momentum dependent, and are given 
in table 1 for the bins of x„. 
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OPAL 




M(7U^7i") (GeV) M(K^K") (GeV) 



Figure 1: Sum of the fitted mass spectra of 7r~^7r~ (with hke-sign spectra subtracted) and K+K" 
for Xp > 0.14 with the results of the fits described in the text. The points show the data, and the 
sohd histograms give the results of the fits. The shaded histograms show the fo contributions, 
the dotted histograms give the f2 and contributions, and the dot-dash histograms show the 
fitted backgrounds. 
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Figure 2: Mass spectra of vr+vr" (with like-sign spectra subtracted) for 6 bins of Xp with the 
results of the fits described in the text. The points show the data, and the solid histograms give 
the results of the fits. The shaded histograms show the fo contributions, the dotted histograms 
show the f2 contributions, and the dot-dash histograms show the fitted backgrounds. 
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OPAL 
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T — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




1 1.05 1.1 1 1.05 1.1 

M(K^K") (GeV) 



Figure 3: Mass spectra of K+K~ for 6 bins of Xp with the results of the fits described in the 
text. The points show the data, and the sohd histograms give the results of the fits. The 
shaded histograms show the fo contributions, the dotted histograms show the (j) contributions, 
and the dot-dash histograms show the fitted backgrounds. 
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Figure 4: Measured fragmentation functions for fg, f2 and (j) together with the output of the 
JETSET 7.4 generator. The f2 and (j) measurements have been scaled by xO.l and xO.Ol 
respectively. The Monte Carlo curves have been normalized to the same total rate as in the 
data. For fo and f2, the errors are statistical only. For the (j), the uncorrelated systematic errors 
are included. There are additional, fully correlated, systematic errors of 7.6% for the fo, 11.6% 
for the f2 and 3.5% for the 0. The x-coordinates for plotting the data points were evaluated 



following |31 . 
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Figure 5: Fitted tt+tt" spectra (with like-sign spectra subtracted) in the fg mass region for bins 
of event charged-particle muhiphcity. The points show the data, and the sohd histograms give 
the results of the fits. The dashed histograms show the fitted backgrounds. 
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Figure 6: Ratios of production rates in data compared to the detector-level JETSET 7.4 model, 
for bins of event charged-particle multiplicity. The errors are statistical only. The ratios have 
been normalized such that the weighted average is unity for each particle. 
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Figure 7: Ratios of total production rates in data compared to the detector-level JETSET 
7.4 model, for bins of the absolute value of the rapidity difference between the meson and the 
nearest charged particle. The errors are statistical only. The ratios have been normalized such 
that the weighted average is unity for each particle. The lowest bin extends down to zero, while 
the highest extends to the kinematic limit. 
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Figure 8: Ratios of total production rates in data compared to the detector-level JETSET 7.4 
model for the energy-ordered jets in three-jet events. The errors are statistical only. The ratios 
have been normalized for each particle to give a weighted average value of unity over the three 
jets. 
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